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Abstract 

We present a description of the freezing of spinodal decomposition in systems, 
which contain simultaneous irreversible chemical reactions, in the hydrody- 
namic limit approximation. From own results we conclude, that the chemical 
reaction leads to an onset of spinodal decomposition also in the case of an 
initial system which is completely miscible and can lead to an extreme retar- 
dation of the dynamics of the spinodal decomposition, with the probability 
of a general freezing of this process, which can be experimetally observed in 

simultaneous IPN formation. 
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While equilibrium thermodynamics describes the overall tendency for a system to stay 
miscible or to realize a phase separation it is dynamics which allows one determine the time 
behavior of a possible phase separation, which depends on microscopic or mesoscopic param- 
eters (e.g. self and interdiffusion coefficients, binary thermodynamic interaction parameters 
and the kinetic coefficients of possible underlying chemical reactions). In general two types 
of phase separation dynamics are expected: spinodal decomposition (SD)! and nucleation 
and growth!. 

Whereas the latter remains effective for the large time regime of the phase separation, the 
initial regime is mainly determined by SD, i.e. nucleation generally precedes the spinodal 
decomposition. Hence, there is no nucleation if a quench of spinodal decomposition is done 
through the critical point. 

A first mean field description of the SD in metallic alloys was given byi, and subsequently 
applied to the SD in binary polymer mixtures by! and!, and is mainly determined by the 
influence of the complicated dynamical properties of the polymer chains on the unmixing 
kinetics. 

In the following we investigate the SD with underlying irreversible chemical reactions, 
which changes both the composition ratio of the compounds and the effective (averaged) 
micro (meso-)scopic parameter of the physical dynamics. Typical cases are noninterfering 
simultaneous polymerization reactions of two monomer types (A and B, called components) 
with a characteristic mutually repulsive interaction, which are involved in interpenetrating 
polymer networks (IPNs)H (or semi- IPNsI) formation. Such growth processes have the char- 
acteristic behavior, that with increasing time the average size of the molecules increases 
(in the case of a percolation-like gelation it becomes infinite for a finite -critical- time) and 
therefore the diffusion coefficients decreases rapidly, i.e. the dynamics of the phase sepa- 
ration becomes very slow and possibly produces an irreversible freezing of the SD in an 
stationary state. Chemical preparation techniques^ and following evaluations using neutron 
scattering! or X-ray scattering^ give at least a qualitative suggestion that the formation of 
homogeneous IPNs can be interpreted as a freezing of spinodal decomposition by irreversible 



chemical reactions, see alsci. Unfortunately, quantitative, time-dependent experiments are 
missing up to now. 

We describe the formation of IPNs (or equivalent systems) by the use of volume fractions 
(pf, where a determines the component (A or B), i characterizes the different molecular 
clusters of one component. In many cases it is sufficient to use as % the number of monomers 
of the different clusters, e.g. tpf is the volume fraction of all clusters with i monomers of 
type A. For each compound (a,i) the mass balance law holds in the hydrodynamic limit 

-J^ + VJ? = Qf (1) 
subject to the conservation of the currents (given by the physical dynamics) X] J° = 0, the 

i,a 

mass of both compounds <pf = f a =constant (f 1 + f 2 = 1) and the source terms Qt — 

i i 

(determined by the chemical reactions under consideration of the noninterfering reactions 
between A and B components). Typical (but not necessarily so) are quadratic forms in the (pf 
for the source terms, which describe reactions of the type [i] + [j] — > [i + j] (cluster-cluster 
react ion0'@) , i.e. Qf = ~ ^fkf'jfk ~ ffY^K^ifJ (for different kinetic coefficients 

j+k=i j 

SO e@0). The current of each compound J° can be described in a linear theory by the well 
known relation J" = — A"iV/if (note, that we use here Einstein's sum convention) with 
the chemical potential fij of the compound (/3,j) and the generalized Onsager coefficients 
(operators) A, which lead in Fourier space to the representation Jf (q) = — iqA?i(q)/^(q). 
Following the arguments ofi, A is defined by A°^ = 5o5{D^S^(q) iff with the diffusion 
coefficient Df and the static structure factor Sf (q) of the compound (a,i), respectively (<pf 
is the averaged (probably time-dependent) volume fraction of this compound). The chemical 
potential follows from the free energy density of a mixture^! 

^ = E[$Wn + f|v<| 2 l+^ £ xT^>l (2) 
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a,(j,n,m 



with the Flory-Huggins parameters which describe the effective binary interaction 

between the components (a,n) and (/3,m). Because this interaction is mainly determined 
by the monomers (and only influenced to a small degree by the molecular topology), we 



can use the more convenient simplification x^s = Xap- gf is a characteristic length, related 
to the subunit length £ of the monomers by ~ ^o/Vn- Hence, using (0), it follows 
immediately that the chemical potential /i° = dF/d(p®- If the deviation of <pf from the 
average iff is sufficient small, we can use the standard approximation of a linearized theory. 
Note that we derive basically a mean field rate law based solely on an assumption that the 
error made in linearization about the instantaneous average is small. As a consequence our 
theory applies only to effects in the lowest order of the kinetic coefficients and the interaction 
parameters Xap- On the other hand, complicated polymeric systems can be described by 
a mean field theory with a sufficient high accuracy^!. However, we get from ([!]) for the 
difference £f = cpf — cpf a Langevin equation 



1 + 2gq 2 



C + E Xa P E C \ + E R nmC + € (3) 
rn 



(g ~ £ ) in which we have supplemented the right hand side by a thermodynamic random 
force 7/" (Rnm follows from the source term by linearization, i.e. R° m = dQ^/dip^ |^=^). The 
correlations of the stochastic forces could be a white noise but this is not necessary as long 
as a spectral function characterizing it exists. The kinetic equation for is the 0-th order 
of the linearization of (jl]) and satisfies d/dt ip° = ~Tp). Thus, (|3]) becomes an internal 
time dependence, because the quantities (p% and R° m are now time dependent. However, 
in the case of an effective binary repulsion between the A and B components, we can deal 
with one effective x-parameter, defined by xab = Xba = Xj whereas xaa = Xbb = can be 
assumed without any restriction. Using the ansatz = (p~%(h a + A"), with the restriction 
J2 tfn^n = 0; i- e - h a = (/ a ) _1 Z)^n (note, that because of the mass conservation the 

n n 

sum over all ip® of one component is a time independent constant), we get from (||) after 
summation over all n: 

h a = {XF (W^(q)> - (G a (q))} h a - (G Q (q)A Q ) + ( V a ) (4) 

with the average over all molecular clusters of one component (...) = 1/ f a ■ ■ ■ an d the 

n 

functions G£(q) =q 2 D«S%(q) [rT 1 + 2gq 2 ] and W£(q) =q 2 D°S°(q). Note, that the source 



term gives a vanishing contribution, because of mass conservation from which it follows 
immediately that S-R° m = dj /dip^J^Qn = 0- The total fluctuation of the composition is 

n n 

given by A£ = J2£n ~ Z)£n = h 1 f 1 — h 2 f 2 , whereas ¥n = 1 an d this leads to h 1 f 1 = 

n n ct,n 

—h 2 f 2 . Hence, 

^ + u{t, q) [T(t, q) - *] A£ = iV(q, *) + I(q, *) (5) 

with u(t, q)- 1 = (2/ 1 (W^q)))" 1 + (2f 2 (W 2 (q))y\ 
N(q, t) = u(t, q) [(r/> / (W l {q)) - ( V 2 ) / (W 2 (q))} and 

I{q,t) = u{t,q) [(G 2 (q)A 2 ) / (W 2 {q)) - (G 1 {q)A 1 ) / (W 1 (q))}. The function 

m i ( Gl (q)) (g 2 (q)) (6] 

{ > q) 2/MW /1 (q)) 2P(W 2 (q)} {) 
can be interpreted as a (time-dependent) kinetic spinodal curve, which describes the onset 
of the SD. Successive approximation, using (H) and (|3|) leads to a perturbation like theory 
for the computation of the (G a (q)A a ). However, the right hand side of (|^) plays only a sec- 
ondary role (especially in the case of symmetric systems (equal chemical kinetic and physical 
parameters) I(q, t) vanishes up to the order of contributions from the thermodynamic noise). 
Thus, (|5]) without the term I(q, t) can used as a preaveraged equation for the SD in systems 
with irreversible chemical reactions. However, the main interest is in the homogeneous part 
of (H), which determines the Greens-function and therefore the evolution of a small fluctua- 
tion. In the case of T(q,t) > x an initially given fluctuation decays exponentially up to the 
order of the thermodynamic noise, whereas for T(q, t) < % an exponential increase of the 
fluctuation occurs, which heralds the SD. Because of the time dependence of T(q,t), we can 
determine the start time t*(q) of the SD (for the mode q) as the solution of T(q,t*(q)) = \ 
(In the case, that for all times T(q, t) < \ is valid, we have t* — 0). Hence, the amplitude of 
the composition fluctuation mode A£(g, t) scales for t > t* as 

|A£(g,t)| 2 ^exp J2 j u(T,q)[x-T(q,T)}dr\ (7) 

This equation allows the determination of the evolution of the composition fluctuations (for a 
given mode q or after integration over all modes q; this is the local composition fluctuations) 



arising as a result of the chemical kinetics and physical (molecular) dynamics. Using the 
usual approximation S%(q) = n[l + |(i?"g) 2 ]~ 1 is the radius of gyration of the molecules 
of the molecular clusters (a,n), which scales as R% ~ n' 7 , the exponent o is determined 
by the underlying model for the structure of the molecular clusters, e.g. in the case of a 
percolation like growth we have a = 0.395 for dimension d — 3) and = DQn" b (e.g. 
b = 1 for the Rouse mo del0, b = 2 for a reptation like mo dell of diffusion), for t — > oo 
it follows u(t,q) — > and T(q, t) — > 2gq 2 ,i.e. the exponent in (0) becomes finite as t — > oo 
(Clearly, this is always the case, if at a finite critical time t c an infinitely large cluster can 
be observed) or at least |A£(g, t)\ 2 no longer scales as an exponential increase exp(ct). The 
importance of this effect is clear: either the SD experiences a dramatic retardation of the 
SD and therefore a later beginning of the subsequent regime of nucleation and growth or 
the SD even freezes into a stationary state without a subsequent segregation regime. This 
can actually be observed e.g. in the simultaneous formation of interpenetrating networksii, 
which leads finally to an averaged domain size in the nm-scale., i.e. the domains can at 
best be interpreted as microphases. This freezing effect can also observed by numerical 
simulations^!. Finally we conclude, that the influence of an irreversible chemical reaction 
on the onset of the evolution of the SD can be sufficiently strong, as seen by: 

1. If in the beginning of the reaction no SD occurs (i.e. for small enough x-parameter) , 
there always exists for sufficiently small q a time t*(q), such that for t > t*(q) the 
SD remains effective. Thus an system which growth by irreversible chemical reactions 
undergoes always the SD. 

2. Furthermore, the same irreversible chemical processes, which induces the SD are in 
competition to the SD, i.e. the existence of the irreversible, simultaneous chemical 
reaction on the thermodynamic separation procedure retards extremely the SD process 
and leads finally to a freezing of the SD, so that irreversible fixed structures are formed. 
A further decomposition after the retardation time scale is controlled by elastic forces, 
which prevent the continuation of the SD. Therefore, the final product contains mostly 
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a strong pronounced elastic tension field. The effect of the chemical reactions on the 
coarsening exponent in ([/]) after the retarded SD starts can be obtained by taking into 
account these increasing elastic forces. Generally, it can be expected, that because 
of these strong forces the contributions of the retarded SD are relatively small in 
comparison to the SD before retardation. 
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